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SPINNING MODE SOUND PROPAGATION IN 


DUCTS WITH ACOUSTIC TREATMENT 
by Edward J. Rice 
Lewis Research Center 

SUMMARY 

Recent acoustic data have shown larger noise attenuations than pre- 
dicted for acoustically treated aircraft engine inlets without splitter 
rings. These encouraging data have stimulated a more detailed theo- 
retical study of the acoustic propagation of spinning modes in acousti- 
cally treated open circular ducts. In addition, the suppressor with 
splitter rings was modeled by using the rectangular approximation to 
the annular duct. The theoretical models were used to determine opti- 
mum impedance and maximum attenuation for several spinning lobe 
numbers from 0 to 50. Some interesting results of the analysis are 
that for circular ducts the maximum possible attenuation and the opti- 
mum wall impedance are strong functions of the lobe number. For an- 
nular ducts the attenuation and optimum wall impedance are insensitive 
to the spinning lobe number for well cut-on modes. The above results 
help explain why suppressors with splitter rings have been quite effec- 
tive in spite of the lack of detailed information on the noise source 
modal structure. Conversely, effective use of outer wall treatment 
alone will require expanded knowledge of the noise source structure. 
Approximate solutions are presented to help interpret the more exact 
theoretical results. 


INTRODUCTION 

Noise suppressor design procedures have evolved over the past 
several years providing quite predictable suppressor configurations 
for aircraft engine inlet and exhaust ducts with splitter rings. These 
procedures have evolved from empirical correlation of extensive noise 
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attenuation tests in flow ducts and from theoretical studies of sound 
propagation in ducts with acoustic treatment. The sound propagation 
theory has used quite simple noise input functions such as a plane wave 
(refs. 1 to 3) or the axisymmetric least-attenuated mode (ref. 4). 

These simple models have provided very good design procedures as 
long as splitter rings were used. Not much attention was paid to the 
suspected, but largely undefined, complex modal structure of the fan 
noise source, and it really did not seem to make much difference that 
this complexity was ignored. It has been known for several years that 
spinning duct modes would be produced by rotor- stator interaction 
(ref. 5) and by inlet distortion- rotor interaction (refs. 6 and 7). Yet 
flow duct tests which did not contain this type of mode produced empiri- 
cal design procedures that resulted in predictable splitter ring suppres- 
sors. 

Recent noise data from aircraft engine inlets showing surprisingly 
high suppression with only outer wall acoustic treatment as well as the 
strong incentive to reduce noise suppressor pressure losses have stim- 
ulated renewed interest in the role of modal structure in noise suppres- 
sor theory. Neither axisymmetric sound propagation theory nor flow 
duct empericism has explained the large noise attenuations which were 
observed. It was hypothesized that the spinning mode character of 
sound does indeed make a difference in the performance of inlet sup- 
pressors with treatment on the outer wall alone. 

The propagation of spinning modes in soft-walled ducts has been 
studied in the past (refs. 8 to 10 are examples), but a general design 
procedure has not resulted from these studies. It appears that a thor- 
ough theoretical and experimental study of spinning mode sound propa- 
gation in acoustically treated circular ducts is needed. This paper 
represents an attempt to initiate this theoretical study to provide rules 
for optimization of duct treatment and to explain some of the observed 
differences between splitter ring designs and outer-wall treatment alone 
configurations. 
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SYMBOLS 

A constant multiplier (see eq. (9)) 
a one- half of lined duct height, m 

optimum resistance coefficient for circular ducts (see eq. (36)) 

B optimum reactance coefficient for circular ducts (see eq. (37)) 

A 

c speed of sound, m/sec 

D circular duct diameter, m 

D m attenuation coefficient (see eq. (60)) 

E acoustic power, W 

f frequency, Hz 

H height between lined walls, m 

J Bessel function of first kind of order m 
m 

k integer variable 

L acoustically treated duct length, m 
M uniform steady flow Mach number 
M t lobe pattern tip Mach number 
m lobe pattern number 

P part of acoustic pressure which is function of transverse coordi- 

o 

nates (r, $ or y, z) (see eq. (2)), N/m 

P part of acoustic pressure which is function of coordinate between 
acoustically treated surfaces (r or y), N/m^ 

2 

p acoustic pressure, N/m 

Q 1 + iM(a + ir) 

R amplitude of eigenvalue a 

r radial coordinate, m 

r^ inner radius of an annulus, m 

r outer radius of an annulus, m 

o 
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t time, sec 

V axial group velocity, m/sec 
gx 

V„ axial phase velocity, m/sec 

px 

v acoustic particle velocity (vector), m/sec 

x axial coordinate, m 

Y Bessel function of second kind, of order m 
m 

y transverse coordinate (see fig. 3), m 

2 

Z liner acoustic impedance (see eq. (11)), kg/m /sec 
z transverse coordinate (see fig. 3), m 

a complex eigenvalue {a = Re 1<p ) 

/3 transverse dimension (either D or H) 

6 radius ratio, r./r Q 

£ dimensionless acoustic impedance, Z/pc 

7 ] frequency parameter, fH/ c or fD/ c 

9 specific acoustic resistance, real (£) 

O 

p density, kg/ m 

cr attenuation coefficient 

r propagation coefficient 

$ angular coordinate, rad 

<p phase of eigenvalue 

X specific acoustic reactance, imag (0 

co circular frequency, rad/sec 

Subscripts: 

j radial mode 

r radial component 

w value at the wall 



5 


INLET SUPPRESSOR DATA COMPARED WITH THEORY 

The motivation for this study came from observation of experimen- 
tal data from inlet suppressors without splitter rings. Samples of these 
data are shown in figures 1 and 2. 

Sound power attenuation spectra from suppressors tested with NASA 
Lewis Fan C are shown in figure 1. Fan C is a high tip speed (472 m/ 
sec, 1550 fps), 1.74-meter (68.3-in. ) diameter, single fan stage with 
26 rotor blades and 60 stator blades (ref. 11). Two attenuation spectra 
are shown: wall treatment plus three splitter rings; and wall treatment 
only. Two interesting points about the data should be noted. First, the 
wall treatment alone configuration is very effective. At high frequen- 
cies, near the blade passage frequency («2000 Hz), the large attenuation 
(~8 dB) obtained cannot be explained by axisymmetric sound propagation 
theory. Second, at low frequencies where multiple pure tones dominate 
(«800 Hz peak), not only is the large attenuation unexplainable by axi- 
symmetric theory, but the splitter rings make very little difference in 
the attenuation. The insignificant effect of the rings for this type of 
sound is probably due to the radial pressure distribution of the sound. 
These questions will be considered again after the spinning mode propa- 
gation theory is developed herein. 

A number of normalized peak attenuation data points are shown in 
figure 2 along with three theoretical curves. The maximum sound 
power attenuation is seen to be normalized by the quantity L/0. The 
length L is the length of acoustic treatment while /3 is a radial dimen- 
sion. The dimension / 3 is the duct diameter for open circular ducts 
and the duct height between treated surfaces for annular or rectangular 
ducts. Note that j8 occurs also in the abscissa (frequency parameter) 
and should be similarly interpreted there. Maximum sound power at- 
tenuation will be explained in more detail later (or see also ref. 1). 

For now, it can be simply considered as the normalized peak experi- 
mental or theoretical sound power attenuation. The theoretical maxi- 
mum attenuation is obtained with the walls having a particular acoustic 
impedance (optimum impedance) at a particular sound frequency. This 
optimum impedance may differ for various duct geometries or for differ- 
ent modal distributions which will be discussed in more detail later. 
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Two of the theoretical curves represent the largest possible attenu- 
ation for circular and rectangular ducts when only the least- attenuated 
mode is present in these ducts. The third curve shows the maximum 
attenuation if a plane wave were impressed upon the liner at one end. 

The least- attenuated mode, often thought of as the most conservative 
indication of liner performance, is seen to potentially damp faster than 
the plane wave input. A combination of modes must be excited to match 
this plane wave input in a soft-walled duct (since a plane wave is not a 
soft wall solution), and this conglomerate is less readily attenuated than 
the so-called least-attenuated mode. 

It should be noted that the curves of figure 2 do not represent an at- 
tenuation spectrum which could be obtained in practice. Each point on 
the curve represents a liner with optimized wall impedance. No real 
liner could track this optimum impedance over a frequency range. At 
most, a real liner could reach these theoretical limits at one frequency 
only. The theoretical curves were generated without grazing flow (M = 0) 
but they will not change with flow in the region where the attenuation is 
falling off with 77 (ref. 3). This lack of dependence on M will be shown 
in a later section. The L/H = 4 plane wave solution was shown since 
it is typical of the splitter ring suppressors tested. Small changes would 
occur for other values of L/H (ref. 3). 

Some interesting observations can be made by comparing the data 
with their respective theoretical limits. First, the rectangular flow 
duct data agree quite well with the rectangular duct theory. The theo- 
retical curve in this frequency range is not influenced very much by the 
type of modal distribution (plane wave, equal energy, equal amplitude, 
etc. ) existing in the duct. This is due to the small number of propagating 
modes which could exist over this frequency range (i. e. , from two to 
five modes). 

The second set of data (engine or fan inlet with splitters) might also 
be compared to the rectangular duct theory. A rectangular duct should 
be a fairly good approximation to at least the outer annulus (which should 
carry most of the acoustic power) of a splitter ring suppressor which 
typically has a radius ratio of about 0.8. These attenuation values are 
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seen to be quite low (about one-fourth) with respect to the theory. 

There were several reasons why this difference between theory and 
experiment may have occurred. The first was intentional; the two sides 
of the annular passage were designed for different frequencies so as to 
broaden the bandwidth of sound attenuation. Thus, at any given frequency 
the entire treated area could not be working at optimum as was the liner 
for the theoretical curves. Several uncontrollable phenomena were prob- 
ably also degrading the liner performance. The liner impedance models 
are not exact, and it would indeed be difficult to construct a liner which 
was truly optimum. Also, the boundary layer refraction effect was not 
considered in the liner design. Whatever the reason, it is not too sur- 
prising to find experimental attenuation levels well below a theoretically 
predicted maximum. It is logical to consider the theory as an ideal up- 
per limit and to use theoretical trends for extrapolation from the ob- 
served experimental data. This procedure has subsequently produced 
quite predictable liner designs where splitter rings are used. 

The third set of data (wall treatment only) is seen to be an order of 
magnitude above the axisymmetric propagation theoretical estimate for 
circular ducts. It is probable that data could be obtained which would 
be even higher if the walls had been truly optimized. These data are 
quite easily explained if one accepts that the sound is mainly concentrated 
in spinning modes, since these modes damp faster than the axisymmetric 
mode theory of figure 2 would suggest. However, spinning modes should 
also exist with splitter rings; and at first one might expect that these 
modes should also damp faster than the axisymmetric modes. The data 
of course deny this conclusion. It will be shown after the development 
of the spinning mode theory that the spinning mode properties are largely 
irrelevant in the splitter ring case, but are the crucial factor in the case 
of wall-treatment only. 

SPINNING MODE PROPAGATION THEORY 


Certain assumptions are common in all of the theoretical develop- 
ment which follows. The steady flow is uniform everywhere in the duct. 
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Boundary layers are sufficiently small so that their effect may be ne- 
glected. The soft walls are assumed point reacting and of uniform im- 
pedance. These assumptions imply that no waves are propagating within 
the liner backing cavity (prevented by honeycomb cells), and that any 
perforations are closely spaced compared to a wavelength of sound. No 
end reflections are considered. That is, the soft duct is sufficiently 
long to render the reflections at the duct termination unimportant; or 
the mode is well enough cut on so that no end reflection will occur. A 
cut- on mode is defined here as a propagating mode; the frequency pa- 
rameter ( 77 ) is sufficiently high so that axial propagation occurs. Only 
very near cut-off will this assumption (no end reflection) be a problem. 
Below cut-off the modes carry very little energy due to their out-of- 
phase pressure-velocity relationship and well above cut-off the end im- 
pedance (ref. 12) is such that very little reflection will occur. It is fur- 
ther assumed that the sound pressure level is sufficiently low so that the 
linearized acoustic wave equation is valid. 

The acoustic wave equation is expressed as: 

v2 p = J. is. + m jk + m 2 A a) 

C 2 at 2 c ® x 3t 3X 2 

The geometric notation used is shown in figure 3. In the following para- 
graphs, the V operator will be defined for Cartesian and cylindrical 
coordinate systems. If separation of variables solutions are sought in 
the form of 


P = 



i cut- (u>/ c )(<7.+iT. )x 
P.e ] 1 


( 2 ) 


where Pj is a function of the transverse dimensions (y, z or r, $) then, 
from the linearized momentum equations, the acoustic particle velocity 
can be expressed as: 
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iVp. 

v. = *- 

J pwQ. 


( 3 ) 


where 


J 


Qj = 1 + iM(Cj + it.) 


(4) 


the j subscript is used to denote the infinite number of possible modal 
solutions. 

Circular Duct 


For a circular duct (cylindrical coordinates), the left side of equa- 
tion (1) can be expressed as: 


V 2 p = i — (r — 

r 3r V drj Jl 


3 2 P + 9 2 P 

a $ 2 ax 2 


(5) 


The remaining term in equation (2) is then sought in the form 


Pj(r, $) =^(r)e 


-im$ 


Thus, using equations (1), (2), (5), and (6) results in 


where 


d 2 ^ dcP- 


J. + r — i + 


dr 


dr 


i'a.r N 

J 


o 


m 


^ = 0 


( 6 ) 


( 7 ) 




00 


[(1 - M 2 )^ + ir.) 2 + i2M(cr. + ir.) + l] 


(8) 


is the equation that links the eigenvalues (a^) with the axial attenuation 
and propagation coefficients (a, r). Equation (7) is recognized as 
Bessel's differential equation with solution 
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/a.r\ a.r\ 

^( r ) = j m - J -i +A i Y -l- J - 




J m 


(9) 


where the oij and Aj must be determined by radial boundary conditions. 
A second constant of integration which usually appears before J m has 
been set equal to 1. 

Equation (9) is valid for annular or open circular ducts. However, 
for added simplicity in this paper, the annular duct case (splitter rings) 
will be approximated by a rectangular duct development in the next sec- 
tion. Of immediate interest, however, is the case of an open circular 
duct. For this geometry, since Y m (Q!jr/r o ) becomes infinite as r — 0, 
the condition must be imposed that 


Aj = 0 (all j) 


( 10 ) 


The remaining radial boundary condition involves the wall impedance at 
r = r Q . The acoustic impedance can be expressed as: 


Z = 



v 


rw 


r=r. 


(ID 


where v rw is the radial velocity within the wall (axial velocity within 
the wall is assumed to be zero). Using continuity of radial particle dis- 
placement the radial particle velocity just outside the wall for the j“ 
radial mode is (ref. 14, p. 314, eq. (34)): 



Q.v 
J rw] 


( 12 ) 


Note that particle velocity continuity is not satisfied unless Qj = 1 
(M = 0). This is due to the artificial slip condition of the fluid at the 
wall. Using equations (2), (3), (6), (11), and (12) and defining a new 
normalized acoustic impedance there results 



( 13 ) 


where 



pc 






c 


(14) 


is the frequency parameter and D is the duct diameter. Using equa- 
tions (9) and (10), equation (13) becomes: 


? = 


i»nQ 1 J m(<V 




(15) 


from which the eigenvalues (a.) can be determined for given wall im- 
pedance, lobe number (m), steady flow Mach number, duct size and 
frequency. If continuity of particle velocity had been used, Qj would 
have appeared to the first power in equation (15). 

The remaining necessary boundary conditions have been subtly 
satisfied due to the choice of solution forms in equations (2) and (6). 
Tangential continuity of pressure and velocity are satisfied if m is an 
integer in equation (6). Only positive x traveling waves are considered 
by the choice of the x function in equation (2); end reflections are ne- 
glected. 


Rectangular Duct 

Due to the relative simplicity of the sinusoidal functions compared 
to the Bessel functions, the annular duct (splitter ring case) will be 
approximated by a semi-infinite rectangular duct. 

For a rectangular coordinate system the left side of equation (1) 
becomes: 
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v 2 p = if + is. + (16) 

Sx 2 Sy 2 3z 2 

The analogy to equation (6) is 

-imz/r 

Pj(y,z) = ^j(y)e (17) 

where z/r Q has taken the place of The spinning wave pattern will 
repeat continuously over z distances equal to the outer circumference 
of the annulus. 

As previously mentioned this rectangular duct solution should be a 
good approximation to the true annular solution if the radius of curva- 
ture is large and if the radius ratio is sufficiently close to 1. The condi- 
tions are usually satisfied at least for the outer passage of a splitter ring 
inlet which passage might be expected to carry most of the acoustic power. 
Using equations (16), (17), and (2) in (1) results in: 


where 



= 0 
J 


(18) 




(1 - M 2 )(<j. + ir.) 2 + i2M(cr. + ir.) + 1 

J J J J 



(19) 


which links the eigenvalue (cr) to the axial attenuation and propagation 
coefficients. Equation (18) has solutions 


c^.(y) = cos 

J 




( 20 ) 


representing symmetric and asymmetric modes, respectively. Follow- 
ing a procedure analogous to equations (11) to (15) provides the wall im- 
pedance equations: 



13 


c = 


i77 7]Qf cos a. 

J J 


a. sin a. 

J J 


( 21 ) 


for symmetric modes, and 


5 = 


-itfTjQf sin ot. 

J r _ J. 

a. cos a. 

J J 


( 22 ) 


for asymmetric modes where the impedance condition was applied at 
y = ±a and where now 


coa _ fH 

77C C 


(23) 


with H being the duct height between the two walls with acoustic treat- 
ment. Again equations (21) or (22) will provide the eigenvalues (Oj) 
when £, m, M, H, and f are provided. 


Acoustic Power Attenuation 

In the more general development where multimodal propagation is 
considered, the acoustic power attenuation is quite complicated (refs. 

1 and 3). The complete pressure and the three velocity components 
must be constructed from the combination of all of the modes. The 
modal participation is determined by the input waveform from the noise 
source. The acoustic intensity is then determined and, since products 
of pressure and velocity are involved, modal cross coupling must be 
considered. The intensity must be integrated over the cross section at 
two axial stations, say x = 0 and L since the radial distribution of 
intensity may change with axial distance. Finally these integrated in- 
tensities which form the acoustic power can be used to determine an 
attenuation from: 

A dB = 10 log 



(24) 
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For the case of a single mode the situation is much simpler. Since 
a mode does not change relative shape as it travels down the duct, no 
integration is necessary. The procedure outlined above reduces to 

/ -2coo\L/c\ coct.L 

AdB=10 1og(e 3 J=-8.68 — L- (25) 

or 

= -17.4 77 7 ) 0 . (26) 

L/D 3 

for circular ducts and 

— ^ = -17. 4 Tirp. (27) 

L/H 3 

for rectangular ducts. Although equations (26) and (27) look the same, 
the 77 terms are defined differently, and the Oj will be different for 
rectangular and circular ducts. 

The damping (a.) and propagation (t.) coefficients can be determined 

J J 

from equations (8) and (19) to be: 


-iM + i-\/ 1 -(1 - M 2 ) 


a j + lT j 


for circular ducts and 


h ) 2 

y?r 77 y 


M 2 



(28) 


(29) 


for rectangular ducts. The last term under the radical in (29) was ob- 
tained from the expression 



15 


M. = !!!° = _?I2L_ (30) 

1 me m(l - 5) 

which defines the tip Mach number of the spinning lobe pattern where 
6 is the equivalent annular duct radius ratio. 

The calculation procedure and the results of the above theory for a 
single mode will now be presented followed by a return to the theory to 
obtain approximate solutions which will help explain the exact solutions. 


CALCULATION PROCEDURE 

The calculation procedure for cylindrical ducts will be outlined 
here. The rectangular duct procedure would follow a similar pattern 
except for the added simplification of having only sine and cosine func- 
tions. The method used to handle the Bessel functions with complex 
arguments will also be outlined. 

Calculation of Maximum Attenuation and Optimum Impedance 

The main complications center upon equation (15) which relates the 
infinite number of eigenvalues (an) to the normalized wall acoustic im- 
pedance (0. If a multimodal analysis were made (ref. 1) the impedance 
could be specified and the eigenvalues determined by iteration, series 
expansion, etc. For a single mode an inverse procedure greatly sim- 
plifies the problem. The eigenvalue is specified and the impedance can 
be immediately calculated. Since a major point of this paper was to 
determine the optimum impedance, this procedure was used. Consider- 
able care was taken to insure that the selected eigenvalue belonged to 
the particular mode of interest. 

For the actual calculations the starting point was one step removed 
from equation (15). Since contours of constant attenuation were desired 
and A dB is proportional to <Tj (eq. (26)), a value of on was first as- 
sumed. A series of values of t. were then chosen which now defines 
the complex eigenvalue (an) through equation (8) or (28) which are ver- 
sions of the same expression. The eigenvalue is then inserted into 
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equation (15) and the complex acoustic impedance is determined. The 
values of M, % and m are held fixed for the calculation. These could 
later be varied to perform a complete parametric study. 

A sample of the calculation results can be seen in figure 4 where 
lines of constant oj and t. are shown on the wall impedance plane. 

Note that both components of impedance have been divided by the fre- 
quency parameter 77. This will make the plots similar for any 77 pro- 
vided M = 0 or that the mode is well cut-on if M * 0. This last point 
will become more apparent in a later discussion. The initial calculation 
used a very small cr and a carefully selected to provide a large 
resistance and near-zero reactance (hard wall eigenvalues well known) 
to insure that the calculation was proceeding in the proper mode. The 
choice of 7\ will determine the mode for a fixed value of Oj. The cal- 
culation then proceeds by increasing the value of Oj until the teardrop- 
shaped constant damping curves degenerate to a point. This point then 
represents the optimum impedance for this mode and the value of cr 
determines the maximum possible sound power attenuation (from eq. 
(26)) for the given value of 17. Actually only one value of 77 need be 
considered for the cases of most interest since the eigenvalue (Oj) at 
the optimum will not change with 77 . This fixed a. can then be used 
to determine the maximum damping at any other 77 by using equations 
(28) and (26). The cases of most interest again are for M = 0 or that 
if M * 0 the mode is cut-on . This really means that Qj in equation 
(15) and defined by equation (4) is not a function of 77 but can be approx- 
imated by: 

Q. » — - — (31) 

J 1 + M 

for well cut~on modes or from equation (4) Qj = 1 for M = 0 at any 
value of 77 . This will be derived in the discussion preceding equa- 
tion (54). 

If a. is pushed above the maximum defined above, one will find 
that the results generate a new set of contours collapsing to a point at 
some lower impedance, these contours representing the next higher 
radial mode. 
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The above procedure provides the maximum damping and optimum 

impedance for the least attenuated mode if the starting values of a. and 

r. are properly chosen, 

J 


Calculation of Complex Bessel Functions 

The Bessel functions with complex argument could be calculated by 
the procedures given in reference 1, However, for the multitude of cal- 
culations required in the previous section it is most efficient to inspect 
one calculation before proceeding with the next. This requires the rapid 
communication which can conveniently be obtained from a desk calcula- 
tor, The desk calculator often cannot handle the large subroutines re- 
quired by the methods of reference 1. Thus the simplified procedures 
which follow were developed so that the Bessel functions themselves 
need not be calculated; only the ratio of two Bessel functions need be 
considered. 

Equation (15) may be rewritten: 


? = 


i nrfe. 


a 


m+1 


J 

(a j) 


m 


J {a.) 

m- y 


(32) 


The ratio J m+ i/ J m can be calculated as a single unit without ever 
considering each part. The Bessel function recursion relationship can 
be written (ref. 13): 


J k-ih ) + J k + i (a j )= ^ J k ( “i > 


(33) 


or 


J k- 1 (Q 9 _ 2k _ J k+l (a j ) 

J k(°V “j J k ( °V 


(34) 
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For any a., a sufficiently large value of k can be chosen so that the 
series representation of the Bessel function is rapidly convergent and 
the last ratio in (34) can be expressed as: 


J k + i ( “j) 


a. 

J 


J k (a.) 2(k + 1) 


1 + 



(k + l)(k + 2) 


(35) 


If the first term of (35) is considered with large k it is seen to be ap- 
proximately the inverse of and quite small compared to the first term 
of (34). In fact the second term of (34) could be neglected to start the 
recursion calculation although for this paper equation (35) was used. 
Equation (34) is repeatedly used with decreasing integers k until 
k = m + 1 and at each step the left hand side provides the inverse of 
the last term of (34) needed for the next step. Finally equation (34) 
provides the inverse of the Bessel function ratio needed in equation (32). 
This procedure provides excellent results with k » 4 [ | . 

RESULTS OF SPINNING MODE PROPAGATION THEORY 

The following sections present some results of the previously dis- 
cussed sound propagation theory. The following discussion pertains to 
the least attenuated mode for a given spinning lobe number. 


Circular Duct 

In figure 5 the maximum possible attenuation for the least attenu- 
ated mode is shown as a function of the frequency parameter r /. Sev- 
eral curves of constant lobe number (m = 0 to 50) are shown. Several 
loci of constant lobe pattern tip Mach numbers (M t = 1 to 4) are also 
shown superimposed. This tip Mach number (Mj. = n?7/m) can be easily 
derived from lobe number, geometry, and frequency considerations. 
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The interesting point to note in figure 5 is that the maximum possible 
attenuation (peak attenuation for an optimized system) is a very strong 
function of the pattern lobe number (m), with larger values of m having 
larger values of attenuation. This may be the explanation for the appar- 
ently high experimental sound attenuations shown in figure 2 for circular 
ducts with wall treatment only. These data were high with respect to the 
axisymmetric theory (m = 0) but not necessarily high with respect to the 
spinning lobe theory. 

It should be cautioned that a comparison of the data from figure 2 
with the theory of figure 5 probably does not determine the spinning lobe 
number because the experimental suppressors were not necessarily op- 
timized. Much higher experimental attenuations might have been at- 
tained with some other liner wall constructions. 

The optimum wall impedance which goes along with the attenuations 

of figure 5 is given by: 


e 


x 


_V_" 

(1 + M) 2 

> 

V 

(1 + M) 2 ^ 


JR 

Tirj 


< 1 


(36) 

(37) 


with the coefficients and given in table I for various values of 
lobe number, and where R is obtained from a = Re 1(p . 

The corresponding eigenvalues at the optimum impedance are given 
in table II. Note that as the lobe number increases the optimum acoustic 
resistance (see B 0 ) decreases and the optimum acoustic reactance (see 
B ) moves toward zero (the liner tuned frequency). This would imply 
that with spinning modes the peak attenuation would occur at a higher 
frequency (for a given liner construction) than would the axisymmetric 
mode. The approximations given in equations (36) and (37) are for well- 
cut-on modes with 


A 

7T 7 ] 


< 1 


(38) 
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where R is the magnitude of the eigenvalue for the mode. The locus of 
R / TiTj = 1 is shown in figure 5. The correction for steady flow Mach 
number will become apparent after a later discussion of approximate 
results. 

It is apparent from the results shown in figure 5, by equations (36), 
(37), and table I that, for optimum design of circular duct acoustic sup- 
pression, much information must be known about the noise source. As 
a minimum at least the lobe number of the dominant mode must be 
known. Preferably the distribution of acoustic power within the radial 
modes should also be known, since radial modes other than the least 
attenuated mode would cause changes to be made in the numerical re- 
sults of figure 5 and table I. The general conclusion that sound attenua- 
tion depends upon lobe number would not be changed for radial modes 
other than the least attenuated mode. 

Rectangular Approximation to Annular Ducts 

The results of the rectangular duct theory, intended to approximate 
the annular duct theory, are shown in figure 6. Again the maximum 
possible attenuation normalized by the duct L/H is shown as a function 
of frequency parameter r\ = fH/c. Several lobe tip Mach numbers and 
lobe numbers are shown. The lobe tip Mach number (M^.), frequency 
parameter (??), and lobe number (m) relationship shown on the figure 
comes from equation (30) with radius ratio 6 = r./r = 0.8. 

In contrast to figure 5 for circular ducts, the rectangular duct 
theory shows that for moderately cut-on modes (Mj, > 1.4) the spinning 
lobe theory collapses rapidly to the axisymmetric propagation theory 
(m = 0). Another way of viewing figure 6 is that for the usual range of 
interest for splitter ring suppressors (!<??< 2), little difference in 
attenuation is observed between 20 lobes and zero lobes. Thus, the at- 
tenuation obtained from splitter ring suppressors is relatively insensi- 
tive to spinning lobe number if the modes are well cut- on. 

The optimum wall acoustic impedance for the least attenuated spin- 
ning mode in a rectangular duct can be given by a single set of equations: 
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6 


0.93 ?7 ^ 
(1 + M) 2 


(39) 


X 


~ -0.74 y / 
(1 + M) 2 


^ < 1 

77 T] 


(40) 


J 

which are applicable for any lobe number (m). These equations again 
represent results for the least- attenuated mode. The optimum eigen- 
value was given by R ~ 2. 364 and <p « 28 degrees. 

Thus, if a splitter ring inlet were being designed for the least- 
attenuated spinning mode and if this mode were well cut-on, the optimum 
liner design would be insensitive to the spinning lobe number as would be 
the attenuation. This interesting conclusion helps to explain the experi- 
mental results for splitter-ring inlets shown in figure 2. Although the 
noise may be distributed in spinning lobe patterns, the axisymmetric 
(m = 0) theoretical attenuation may be close to an upper limit for this 
type of suppressor. 

The insensitivity of splitter- ring suppressor design and attenuation 
to spinning lobe number may also explain the quite good results obtained 
from previous suppressors in spite of the meager information available 
about the noise source. 


Comments on Rectangular Approximation to Annular Duct 

The use of the rectangular wave equation (eqs. (1) and (16)) to re- 
place the more exact cylindrical coordinate system of the annular split- 
ter ring passage should create no problem for the limited scope solu- 
tions of interest in this paper. For the high radius ratios of interest, 
curvature effects should be negligible. In reference 14 radius ratios 
from zero to about 0.8 were used in a cylindrical system with constant 
frequency parameter ( 77 ), and no change was found in the propagation or 
attenuation coefficients when r./r Q > 0.4. This implies that curvature 
effects were negligible above this radius ratio (0.4). In reference 14 
only non spinning mode results were shown, and it is conceivable that 
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spinning modes with high spiral angles may be more sensitive to curva- 
ture. However, at least for well- cut-on spinning modes, there is very 
little spiral angle in a cylindrical duct which can be seen by looking at 
the axial phase velocity: 


v « c(l + M) 


or the axial group velocity: 


v gx ^ + 


1 + cos 2„ 

2 \J”)/ 


(41) 


j (1 +M) 
2 


(“f 

w 


cos 2<p\ 


(42) 


which are respectively slightly greater and less than those of a plane 
wave (no spiral angle) in a hard duct which are both equal to |c(l + M)]. 
This is true since (R./tt r?) is small for a well cut- on mode. 

Another problem does arise, however, with regard to what modal 
structure of sound is excited within the annular liner passages . The 
results of this paper are for the least attenuated mode in either the cir- 
cular or annular duct. This represents the most conservative single 
mode result. The problem in the splitter ring passage can be illustrated 
in figure 7. The first radial mode (cylindrical coordinates) with small 
and large lobe numbers are sketched as input to the splitter ring pas- 
sage. For small enough lobe number, the pressure is more uniform 
over the annulus and the least attenuated mode will dominate in the 
carrying of sound power. The results of this paper should be valid for 
this case. As lobe number is increased, more higher radial mode con- 
tent will be excited in the annulus, and the least attenuated mode approx- 
imation will underestimate the true attenuation because these higher 
radial modes have higher damping coefficients. 

If the lobe number is increased even further (dashed line, fig. 7), 
the pressure amplitude at the splitter ring will approach zero, and the 
ring will be completely ignored in the sound propagation. For this case 
the open circular duct (no splitter rings) solution will become valid. 

This can be further emphasized by repeating equation (9) 
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( 9 ) 


which represents an annular duct solution, and observing the value of 
A 1 (first radial mode) in figure 8. Figure 8 results are from refer- 
ence 15 for hard wall annular ducts. For a given radius ratio (5), as 
the lobe number (m) is increased, A 1 rapidly drops toward zero. Of 
course as A.^ - 0, the open circular duct solution is approached, that 
is: 


£J(r) - J 


m 



(43) 


and the presence of the splitter ring becomes irrelevant to the solution. 
The same phenomenon occurs for higher radial modes but at higher lobe 
numbers. This type of modal pressure distribution may account for the 
noise suppressor data at low frequencies shown in figure 1. Recall that 
for frequencies of about 1000 hertz and lower, the splitter rings pro- 
duced very little added attenuation. 


INTERPRETATION OF THEORETICAL RESULTS 
USING APPROXIMATE EXPRESSIONS 

It is useful to derive approximate expressions from the preceding 
theoretical analysis to clarify the basic difference between noise sup- 
pressors with splitter rings and with wall treatment alone. Approxima- 
tions stated earlier without proof will be derived here. The approxima- 
tions are mainly for well cut-on modes, which are of the most interest, 
since they carry the bulk of the acoustic power. 


Well Cut-on Modes 

Some preliminary approximations are required; the first being 
applied to equation (29) for the damping and propagation coefficients 
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for rectangular ducts. The subscripts will be dropped and it will be 
understood that the least attenuated mode will be considered. 



For a well cut-on mode: 

(1 -M 2 ) 


+ - 
\ vr] ) M‘ 


“tj 


« i 


which allows the approximate expression 




x ~ 1 - 


(x « 1) 


to be used. Also use: 

Q! 2 = R 2 e* 2<iP = R 2 (cos 2<p + i sin 2 <p) 
Therefore (29) can be written: 


1 (* y sl , 

2 W 


cr + ir « — ( — sin 2cp + i 


1 1 


II + M 2 


cos 2(p + -L 
M? 


tj 


The damping coefficient is thus: 


„ •my* 

2 \7T?7/ 


sin 2 <p 


(29) 


(44) 


(45) 


(46) 


(47) 


(48) 


and the propagation coefficient is: 
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1 + M 2 



cos 2 (p + 


< 


(49) 


For circular ducts equation (28) can be seen to be the same as equa- 
tion (29) except that the term is missing, and the lobe pattern tip 
Mach number information is in the (ot/irv) 2 term. For circular ducts 
the equations can be immediately determined from (48) and (49) by ne- 
glecting the Mj. term. 



a « — 


sin 2 (p 


(50) 


2 

cos 2 tp (51) 

Note that the Mach number (M) has cancelled out of the damping equa- 
tions (48) and (50). It will be shown shortly that the eigenvalue (« = Re lt() ) 
is constant at the optimum point (limiting contour in fig. 4) even with a 
steady-flow Mach number. Thus, the maximum possible attenuation is 
not a function of Mach number but depends only upon the frequency pa- 
rameter ? 7 . 

Another approximation required is for the Q. term in the wall im- 
pedance equations. This was given in equation (4) as: (again dropping 
subscripts) 


1 + M 



Q = 1 + iM(a + ir) (52) 

Note that in equation (47) the eigenvalue had to be carried to second 
power to get the damping coefficient a. However, for well cut-on modes 
equation (47) is dominated by: 

a + ir w — - — (53) 

1 + M 


thus equation (52) can be estimated by 
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(54) 


Q 


1 + M 



Equation (54) can now be used in the wall impedance equations (21) 
and (15) to provide 

g(l + M) 2 ^ i?r cos a ( 55 ) 

ri a sin a 


for rectangular ducts and 


m + m) 2 

V 




aJ m+l {a) - mJ m (a) 


(56) 


for circular ducts. 

Rectangular approximation to annular ducts . - Some interesting 
conclusions can be drawn from equation (55) for rectangular ducts. 
Consider the left side as an effective impedance and note that the right 
side is a function only of the eigenvalue. Thus, an optimization pro- 
cedure as explained in the CALCULATION PROCEDURE section and as 
shown in figure 4 need be done only once if the coordinates in figure 4 
are considered as effective impedance instead of as impedance only. 

The eigenvalue (a) at which the constant a contours collapse to a point 
(optimum effective impedance) will always be the same. Since a = Re 1<f> 
is always constant at the optimum effective impedance, the conclusion 
stated after equation (51) is now justified. Since the damping coefficient 
(a in eq. (48)) is not a function of Mach number (M) and a is constant, 
the optimum damping is a function only of the frequency parameter tj. 

To carry this conclusion one step further, from equation (27) the nor- 
malized maximum possible attenuation (A dB/L/H) is a function only of 
the frequency parameter rj for well cut-on modes and is given by 

A dB ^ 12.8 

L/tt V 

where the optimum rectangular eigenvalue (R 
used. 


( 57 ) 

= 2. 364, <p = 28°) was 
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Going back to equation (55) if the optimization procedure is per- 
formed at M = 0 and r\ = 1 it is found that 

c = 0.93 -i 0.74 (58) 

thus in general 

g ~ (0.93 - i 0 . 74)77 ( 59 ) 

(1 + M ) 2 

which justifies equations (39) and (40). 

Recall that the lobe number (m) did not occur either in the damping 
coefficient (eq. (48)) or in the wall impedance (eq. (55)). Thus for well 
cut-on modes the maximum possible attenuation and the optimum wall 
impedance are not functions of the spinning lobe number, a result which 
was also shown by the more exact theoretical results. 

Circular ducts. - Some important differences between circular and 
rectangular ducts can be deduced from the approximate equations. For 
circular ducts the right hand side of the effective impedance Eq.' (56) is 
Seen to be a function of the lobe number (m). Thus although the optimiza- 
tion procedure of Fig. 4 can be performed for M = 0 and n = 1 and then ap- 
plied to all M and ^ the procedure must be repeated for each value, 
of lobe number (m) of interest. This results in the optimum impedance 
coefficients of table I and the optimum eigenvalues of table II being func- 
tions of lobe number. This also implies through equations (50) and (26) 
that the normalized maximum possible attenuation (A dB/L/D) is a func- 
tion of lobe number. It is still not a function of steady flow Mach number 
for R /7717 < 1 in figure 5. Equations (27) and (50) can be combined to 
yield 

A dB « (60) 
L/D v 

where D m is an attenuation coefficient dependent upon lobe number (m) 
and is given in table III. 
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Note that equation (60) used in conjunction with table III yields an 
approximation for large values of the frequency parameter ( rj ). This 
inverse function is valid for the portions of the curves in figure 5 which 
are straight lines. 


Cut-off Modes 

When a mode is cut-off it is of less interest since it carries no 
power and is attenuated in even a hard walled duct. For completeness 
of the approximate equations, the attenuation of the cut-off modes will 
be developed. The term in equation (28) or (29) will dominate 

the radical and with M = 0 the result will be: 

a w R cos a (61) 

ITT) 

for either rectangular or circular ducts. 

Using equations (26) or (27) yields 

— « -17.4 R cos a (62) 

L//3 

where |3 = D for circular ducts and 0 = H for rectangular or annular 
ducts. Note that this expression does not contain the frequency param- 
eter 7] and represents the horizontal portions of the curves in figures 
5 and 6. This behavior has been observed before in references 16, 1, 
and 2. For M = 0 the same eigenvalues given following equation (40) 
for rectangular ducts and in table II for circular ducts can be used. 
However for M * 0, the situation is much more complicated. The sim- 
ple expression (eq. (54)) for Q in the wall impedance equation is no 
longer valid. Q will be a function of % and instead of just a simple 
stretching or compression of the equal damping curves in the impedance 
plane (see fig. 4) a rotation will also occur. The eigenvalues will thus 
be changed from those of a well cut-on mode. 

An example of the effect of Mach number on the maximum possible 
attenuation for the least attenuated mode in rectangular ducts can be 
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found in figure III-13 (p. 87) of reference 4. It should be noted that only 
at low values of frequency parameter is this maximum attenuation a 
function of Mach number where attenuation increases for inlets and de- 
creases for exhaust ducts. For 77 > R/tt this Mach number dependence 
vanishes as previously discussed. 

CONCLUDING REMARKS 

This paper could be briefly summarized as follows. Suppressors 
with outer wall treatment alone have been observed to provide much 
more attenuation than previously expected. Such expectations were 
based upon the extrapolation from splitter ring suppressor measure- 
ments or least attenuated axisymmetric, circular duct theory. The 
theoretical analysis of the least attenuated spinning mode with sound ab- 
sorbing walls has revealed that the experimental results could have been 
anticipated qualitatively. Approximate expressions for the attenuation 
of well cut-on modes were presented to illustrate the theoretical results. 
To further pursue this approach more information is required about the 
spinning modal structure of the sound emanating from the source. 

It was consistent with the spirit of this paper, then, to make all 
necessary assumptions required to point out as simply as possible the 
gross behavior of both splitter ring and wall only suppressors. Many 
of the conclusions made here could possibly be challenged by postulating 
modes near cut-off or more complex radial mode structures. However, 
the experimental data do provide credibility, if not precision, to the as- 
sumptions made here. 


CONCLUSIONS 

Within the restrictions of the assumptions made within this paper 
the following conclusions were made: 

1 . Noise suppressors with splitter rings are relatively insensitive 
to the spinning mode lobe number. This explains why quite simple non- 
spinning mode theory and rectangular flow duct tests have provided quite 
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predictable splitter ring designs in the past even where spinning modes 
were possible. 

2. Circular noise suppressors with wall treatment only are very 
sensitive to spinning mode lobe number. This explains why axisymmet- 
ric cylindrical duct theory does not agree well with experimental data 
taken with fan noise sources. 

3. Due to conclusion 2 above, much more detailed information about 
the fan noise source will be required to best take advantage of the great 
potential of wall- treatment-only suppressors. 

4. Approximate expressions for sound attenuation for well cut-on 
modes have been presented which provide insight into the difference be- 
tween the two types of suppressors. 

5. In a splitter ring suppressor (or an annular exhaust duct), the 
presence of the inner wall can under some circumstances be completely 
ignored by the propagating spinning mode. This will occur when the 
radial pressure distribution is skewed toward the outer wall and very 
little pressure oscillation occurs on the inner wall. A high spinning 
lobe number and low radial mode will cause this to occur. 
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TABLE I. - OPTIMUM RESISTANCE AND REACTANCE 
COEFFICIENTS FOR SPINNING MODES 
IN CIRCULAR DUCTS 


m 

0 

1 

— 

3 

5 

7 

10 

20 

50 

B e 

0.890 

0.644 

0.448 

0.357 

0.303 

0.251 

0. 170 

0.097 

B x 

0.38 

0.221 

0. 121 

0.095 

0.076 

0.060 

0.038 

0.020 


TABLE II. - EIGENVALUES AT OPTIMUM IMPEDANCE 
FOR SPINNING MODES IN CIRCULAR DUCTS 




1 

3 

5 

7 

10 

20 

50 

R 

3.278 

4.691 

7.299 

9.759 

12.088 

15.547 

26. 662 

58.620 

<Pt 

deg 

23.19 

18.07 

13.49 

10.91 

9.75 

8.03 

5.46 

3.32 


TABLE m. - COEFFICIENT FOR MAXIMUM POSSIBLE 
ATTENUATION, , SPINNING MODES, CIRCULAR DUCT 


m 

0 

1 

3 

5 

7 

10 

20 

50 

D 

m 

21.51 

35.88 

66.83 

97.87 

134.9 

184.9 

372.3 

1099 
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Figure 1. - Effect of inlet suppressor configuration on sound power 
attenuation, {Fan C, 90% speed, ) 
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Figure 2. - Comparison of suppressor theory and data. 






